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Single three dimensional polymers confined to a slab, i.e. to the region between two parallel plane 
walls, are studied by Monte Carlo simulations. They are described by A''-step walks on a simple 
cubic lattice confined to the region 1 < z < D. The simulations cover both regions D « Rf 
and D » Rf (where Rf ~ A'^'^ is the Flory radius, with u ~ 0.587), as well as the cross-over 
region in between. Chain lengths are up to A'' = 80, 000, slab widths up to D = 120. In order to 
test the analysis program and to check for finite size corrections, we actually studied three different 
models: (a) Ordinary random walks (mimicking 0-polymers); (b) Self-avoiding walks (SAW); and 
(c) Domb- Joyce walks with the self-repulsion tuned to the point where finite size corrections for free 
(unrestricted) chains are minimal. For the simulations we employ the pruned-enriched-Rosenbluth 
method (PERM) with Markovian anticipation. In addition to the partition sum (which gives us a 
direct estimate of the forces exerted onto the walls), we measure the density profiles of monomers 
and of end points transverse to the slab, and the radial extent of the chain parallel to the walls. All 
scaling laws and some of the universal amplitude ratios are compared to theoretical predictions. 



I. INTRODUCTION 

Although the behaviour of flexible polymers in a good 
solvent confined to different geometries has been studied 
for many years P, Q, there are still a number of open 
questions. In the present work we shall only discuss single 
polymer chains between two parallel walls which act only 
as geometric constraints, without any energetic effects. 

Theoretically, this problem is rather well understood. 
All important scaling laws have been formulated, includ- 
ing the cross-over from the region of narrow slabs (where 
the distance D between the walls is smaller than the 
Flory diameter of a free coil) to the opposite case of wide 
slabs. In particular, there exists an important theoretical 
prediction: Near such a wall the monomer density profile 
increases as 

p{z) ^ z^l'^ , (1) 

where z is the distance from the wall and v is the Flory 
exponent [J] . This is supposed to hold for all dimensions 
(not only for c? = 3), and both for ordinary random walks 
(ideal polymers) for which v —Xjl and for self avoiding 
walks with v w 0.587 (in d = 3). 

It is intuitively obvious that the force exerted by the 
polymer onto the wall is proportional to the monomer 
density near the wall. The ratio between the two can be 
expressed in terms of a universal amplitude ratio which 
is easy to calculate for ideal chains, and which was cal- 
culated by Eisenriegler Q as an expansion in e = 4 — d. 
Several authors have tried to verify these detailed pre- 
dictions by Monte Carlo simulations 0, IE S 1 but the 
results are not yet convincing. While the scaling of the 
density near the wall is roughly verified, the amplitude 
ratio consistently has come out too large, casting even 
doubt on the validity of the e-expansion. 

In Ref. we had studied confined polymers in a strip 
in two dimensions where the amplitude ratio had been 
predicted by Eisenriegler (using conformal invariance 



results of Cardy et al. ^3)- There we verified all predic- 
tions, but we found that this was less easy than antici- 
pated: There are very large corrections to Eq. which 
can easily be missed, and overlooking them would give 
wrong estimates of the amplitude ratio. This suggests of 
course that the same effect was the source of difficulties 
in 0? 3. 

It is the purpose of the present paper to present a care- 
ful numerical study, in order to settle these questions. 
We not only simulated much larger systems than previ- 
ous authors, going to chain lengths up to TV = 80, 000, 
slab widths up to D = 120, and collecting rather high 
statistics. Since it is well known that asymptotic scaling 
of unconstrained self-avoiding walks (SAWs) is reached 
rather slowly, with correction terms decreasing only as 
El El m, we studied also the Domb- Joyce (DJ) 
model |l4| with w = 0.6 (where convergence to asymp- 
totia is much faster 0,0|) in addition to SAWs. 

The DJ model is defined by the partition sum 

Zn{w)= (2) 

con figs. 

where the sum extends over all random walk (RW) con- 
figurations with A'' steps, < u> < 1, and k is the to- 
tal number of monomer pairs occupying a common site. 
For w = 1 the DJ model describes just ordinary random 
walks. For w = it is just the SAW model. Asymptoti- 
cally (for ^ oo) the model is in the SAW universality 
class for all w < 1, but the speed with which the renor- 
malization group fixed point is approached depends on 
w. Moreover, it is approached from opposite sides when 
w < w* and when w > w* , with w* « 0.6 E3l| . 
For w = w* the approach to asymptotia is fastest, and 
we thus expect also smaller finite size corrections for the 
present problem of confined polymers. 

Finally, we also performed simulations of the ordinary 
RW model, just to check the simulation and analysis pro- 
grams, as everything can be calculated there analytically. 



2 



For the simulation we used the pruned-enriched- 
Rosenbluth method (PERM) [l^ with fc-step Markovian 
anticipation 0, 0, 0, 0| . Apart from being fast (no- 
tice that the pivot algorithm which is very fast for un- 
constrained polymers is very inefficient for narrow 
slabs), it has the advantage that the partition sum is 
computed by default with very high precision. Thus we 
could estimate the dependence of the monomer fugacity 
on the width D, and from that directly the total forces 
exerted onto the walls. For more details see Ref. . 

Details of the scaling predictions are discussed in Sec. 
II, while results and their comparison with theoretical 
predictions are presented in Sec. III. Conclusions are fi- 
nally given in Sec. IV. 



II. SCALING PREDICTIONS 

The end-to-end distance Rn oi a, free SAW in infinite 
d— dimensional volume scales as 



R 



N 



((xjv - xo)^) « d [kNY^- (1 + b/N'^-) (3) 



where Vd is the Flory exponent, is the leading correc- 
tion to scaling exponent, and h and k are non- universal 
constants which depend on the microscopic realization. 
In d = 2 one has 1^2 = 3/4 H, while the best pub- 
lished estimates for d = 3 are — 0.5877(6) [n|, 
0.5874(2) [3, and 0.58758(7) O. In the latter paper 
it was assumed that the leading correction to Eq. (0) is 
1 + b/N^'^ . Since this is ques tionable (our own simula- 
tions gave A3 K, 0.45 [l^ l2f]l|'l. we use in the following 
1/3 = 0.58765(20) which also incorporates results from 
extensive simulations of the DJ model with w — 0.6 [20| . 
For the DJ model w* is defined by b(w — w*) — 0. The 
absence of large corrections to scaling leads to a rather 
precise estimate of the constant k in case of the D J model 
with w = 0.6: k = 0.3259(4) [13, but for 3-d SAWs the 
estimate is much less stable. The value k — .4640(4) 
of Ref. 10 depends crucially on the estimate A « 0.56 
made by these authors. Assuming instead A = 0.5 as in 
Ref. [ij and the value of 1^3 found by these authors, the 
same data would give k « 0.4655, while an even larger 
value would be obtained if A < 0.5. In the following we 
shall use k = 0.4657(7) for 3-d SAWs. 

Equation Q with d — 2 also describes the behaviour of 
the parallel components (i.e. parallel to the wall) in the 
regime 1 <C -D <C N"^^ where the polymer is essentially 
two-dimensional. The constant k depends then on the 
slab width D. A scaling ansatz for the cross-over between 
the two regimes 1 < N"'' < D and K D < N"'' is [l| 



R%^^^{D)=R% ^Rn/D) 
where Rj^ is given by Eq. Q with d = 3, and 



$(7?) 
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„2(l/2/l'3-l) 



for rj 
for 77 



0, 
00 . 



(4) 



(5) 



It leads to the prediction [g| 



(6) 



The partition sum of a free SAW in infinite volume 
scales for N ^ 00 as 



Mo 



-N 



const 



(7) 



with floe being the critical fugacity per monomer, and 
with being a universal exponent. In two dimen- 
sions 72 — 43/32 01, while the best published esti- 
mate for d = 3 is 73 = 1.1575(6) In the fol- 
lowing we shall use the estimate 73 = 1.1575(3) ob- 
tained from the DJ model with w = 0.6 For 3-d 
SAWs one has fioo = 0.213491(4) from exact enumer- 
ations [23 and/Xgo = 0.2134910(3) from Monte Carlo 
simulations [13.1201 . In the following we shall use the lat- 
ter. For the DJ model we use /Xoo = 0.18812145(7) 

Again we must expect that the same ansatz, with 73 
replaced by 72, with /ioo replaced by fi{D), and with the 
constant replaced by c{D), holds for slabs in the limit 
D < N"^: 

Zn{D) = fi{D)-^N^^-^ c{D) for D < N"'' . (8) 
The corresponding cross-over ansatz is then 

Zn{D) = Zn ^(Rn/D) (9) 

with 



r const for rj ^ , 

*W = I ^(72-73)/i'3 exp(-^(7//^/3)l/''^) for r] ^ 00 . 

(10) 



and 



i{D) = fioo+aD-^^''' 



c{D) ^ const D^^^-^^^/''^ 



(11) 



(12) 



Notice in particular that the Z?-dependence of IJ.{D), 
Eq. Hll|l . follows directly from the scaling ansatz for the 
cross-over. We should also point out that the partition 
function for a polymer in a slab is defined as the sum over 
all walks starting at fixed Xo.n = (xq^uq), but averaged 
over all zq & [1, 

The force exerted onto the wall is most straightfor- 
wardly expressed in terms of the work done when moving 
one of the walls, i.e. by the dependence of the free energy 
- and thus also of the partition sum - on D, 



F = kRT 



d\nZN{D) 
dD 



(13) 



where we have introduced a dummy temperature T which 
can take any positive value. From Eqs. JS)) and Hll|l . the 
force per monomer is then obtained as 



f^F/N^ keT 



£)-l-l/l'3 



(14) 
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in the limit of D ^ oo and N » D'^/"^ . 

The monomer density near the wall is predicted to 
scale with Eq. Q . For ordinary random walks this gives 
p{z) ~ z^, but in that case one can compute p{z) exactly, 
with the result 



(a) 



Dp{z) = 1 - cos(27rz/i:>) = /(z/D) 



(15) 



for D ^ 1 (we normalize p{z) such that X^^i Pi-^) = !)• 
For walks with excluded volume p{z) is not known, but 
one expects Dp{z) to be universal. 

One should expect that the density near the walls is 
proportional to the force per monomer. Indeed it was 
shown by Eisenriegler Q that 



lim k ^ 



= B 



f 



= B- 



D 



-1-1/1/3 



(16) 



with B being a universal amplitude ratio. For ideal 
chains one has i? = 2, while for chains with excluded 
volume in 4 — e dimensions one has i? « 2(1 — 6ie) with 
61 = 0.075 "a^. In three dimensions this gives the predic- 
tion B « 1.85. 

Finally, the density profile of end points scales near 
the walls with a new exponent which is related to the 
surface exponent in spin systems j23j . For polymers, the 
standard way to introduce this exponent is via the par- 
tition sum of a SAW, one end of which is glued to an 
impenetrable wall. For this system one has 



7(1) 



(17) 



with 7^^^ = 0.679(2) |2j|. The end point density then 
scales like [2^ 

p^U^) ^ z(7-7<^>)/. ^ ^0.814(6) (18) 

where the right hand side holds for d = 3. 

III. RESULTS 
A. Random walks 

In order to test our simulation and analysis methods, 
we first simulate the simple model of RWs on a simple 
cubic lattice between two hard walls at z = and z = 
D + 1, confining the polymer to 1 < z < D. We simulated 
width D up to 80 and chain length between 3500 (for 
£) = 4) and 80000 (for £> = 80). 

Monomer densities are shown in Fig. ^ They were 
obtained by averaging over the central part of the chains, 
excluding 10 % on either side to avoid errors from the fact 
that Eq. should hold only far away from the ends, 
for monomer indices n satisfying ^ n ^ N — 
(we should mention that N / D"^ > 10 for all data sets). 
For finite D the scaling has to be slightly modified, by 
replacing in Eq. ifTCjl z/D by ^ = z/{D + 1) and Dp{z) 
by {D + l)p{z). We see that all data, even for small D, 
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FIG. 1: Rescaled values of the monomer density, {D + 1) p{z) 

against ^ = z/{D+l) for ordinary random walks. Also plotted 

is the function /(^) = 2sin^(7r^). (b) The same values, but 
divided by /o(C) = 2^2[^(l^^)j2^ 



fall precisely onto the predicted curve. To show that also 
the regions near the walls are correctly sampled, we plot 
in panel (b) of Fig. ^ the same data but divided by the 
product of the two power laws for z and z k D, 
{D + l)p(z)//o(0 with /o(0 = 2^2[^(i _ ^)]2. 

Critical fugacities are determined by plotting Z^p^ 
against log TV and demanding that these curves become 
horizontal for large N. Results are shown in Fig. [3 where 
we plot pd — Moo with p^a = 1/6. The dashed line is not 
a fit to the data, but fits their extrapolation to Z? — > oo, 
pjj - p^, = 0.2741/(D + 1)2. This agrees with Eq. ^ 
(since — 1/2) and gives a = 0.2741(2), where the error 
is obtained by assuming that the slope is 2 as predicted. 
On the other hand, Eqs. (|16|l and (|15|l together with B ~ 
2, k = 1/3, and = 1/6 give a = 7rV36 = 0.27416, 
in perfect agreement. We should point out that Eq. Hll|) 
is significantly violated for small D in this model, and 
becomes exact only for large D. We shall see the same 
behaviour also for SAWs and for the DJ model, and we 
had seen the same also in d = 2 . 

In Fig. 121 we show the transverse distribution of chain 
ends. For ordinary RWs it should be just the square 
root of the monomer density, i.e. proportional to sin(7r^). 
This is obviously true for our data, with high precision. 

Finally we should mention that the cross-over ansatzes. 
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FIG. 2: Log-log plot of /id — fJ.oo against D + for ordinary 
random walks. The dashed line is fin — Moo = 0.2741/(D + 1)^ 
and gives our best "extrapolation" of the data for large D. 
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FIG. 4: Data collapse for testing the cross-over ansatz Eq. Q 
for self avoiding walks. 
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FIG. 3: Rescaled values of the probability Pend{z) that the 
chain end is at the distance z from a wall against ^ = z/{D + 
1), for ordinary random walks. The solid line is the function 
f sin(C). 



Eqs. Q and Q become trivial for RWs, since the crit- 
ical exponents are the same in c? = 2 and d — 3. We 
therefore also do not show data for -Rjv.|| which can be 
calculated approximately by assuming that all steps are 
uncorrelated, and that vertical steps occur with proba- 
bility 1/3 inside the slab and with probability 1/5 at the 
boundaries. 



B. Self-avoiding walks 

Let us first discuss the trivial case D = 1. In this case 
we have ordinary 2-d SAWs, and therefore we can use the 
comparison with the known results as a test for our algo- 
rithm. Our simulations, with N = 3000, gave indeed 
perfect agreement for the critical exponents, and also 
the amplitude for the end-to-end distance, R^^/N'^'^^ « 
0.771(1), in agreement with the value obtained in [lH . 

In our non-trivial simulations we used widths up to 
D = 120 and chain length up to iV = 80, 000. As first 




0.1 1 10 

Rn/(D+1) 



FIG. 5: Data collapse for testing the cross-over ansatz Eq. @ 
for self avoiding walks. 



tests we checked the cross-over ansatzes Eqs. Q and (O- 
In these tests we replaced Rn and by parametriza- 
tions similar to Eqs. (@J) and @, but including additional 
correction to scaling terms. As seen from Figs. 0] and |S1 
the data collapse is excellent. Such a perfect collapse 
would not have been obtained, if we had replaced i?jv 
and by the leading asymptotic powers of N p^ . 

Critical fugacities were determined by plotting 
log^Af — (72 — l)logiV -I- Nx against N and changing 
X until these curves become horizontal for large N . Then 
^.{D) = exp(a;). Results are shown in Fig. El where we 
plot ijl{D) — jjLrx, against D. As for ordinary RWs, the plot 
does not give a straight line (replacing Dhy D + 1 would 
improve the situation a bit, but not much), so the straight 
line shown in Fig. |S1 indicates the estimated asymptotic 
behaviour, assuming its slope to be given by — l/i/3. It 
provides us with the estimate a — 0.448 ± 0.005. 

The monomer densities for different values of D, again 
from the central region only, are shown in Fig. [7| Plot- 
ting the densities directly as in panel (a) indicates that 
scaling is satisfied. But it is not very informative, since 
deviations from scaling in the important regions near the 
walls would not show up. Also, panel (a) might suggest 
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FIG. 6: Log- log plot oi fio ~ j-ioa against D + 1. The dashed 
line is fio - fiao = 0.448(D + l)"^''''^ 



that p{z) is simply the product of two powers 
1 



p{z) 



r/o( 



D + 1 'D + 1' 



with Mo^Am-o]'^"' , 

(19) 

where the constant A = 18.74 is determined by nor- 
malization. We had already seen for random walks 
that p{z) is not that simple, and indeed plotting {D + 
1) p{z) / fo{z / {D + 1)) as in panels (b) and (c) shows that 
this would be a very bad approximation. In addition, 
panel (b) shows that one has to introduce an 'extrapola- 
tion length' S as suggested in Ref. @, so that the scahng 
variable ^ is replaced by 



6 = 



z + S 



D + 1 + 26 



(20) 



Best scaling near z — and z — D (panel (b)) and best 
data collapse (panel (c)) is obtained for 6 « 0.15, al- 
though a closer inspection of these figures shows that nei- 
ther the scaling nor the data collapse are perfect. These 
small persistent discrepancies and the overestimation of 
the amplitude B discussed in the next paragraph were 
the main reason for studying the Domb- Joyce model. 

FigureEt suggests that D^+^Z''-' p{z)/ z'^/"^ 0.87(4) x 
A = 16.1(8) for z and D oo. The very large uncer- 
tainty reflects the rather steep slopes at z = and z = D. 
Using this in Eq. (0 gives B = 2.13±0.11. This is larger 
than the prediction of Eisenriegler , but much less so 
than previous estimates |^ Q ■ We believe that these pre- 
vious authors had missed the fact that p{z)/ fo{£,) is not 
constant. If we would assume p{z) cx /o(C)i we would 
obtain B « 2.48 which is indeed similar to the previous 
Monte Carlo estimates. 

As a further test of scaling we checked in detail that 
i?ll(D) ^ N"^ for N"^ > D, and we estimated the 
asymptotic ratios between the two. They are plotted 
in Fig. where we also plotted the scaling prediction. 



2 1^2 



(21) 



Finally, we show in Fig. |51the distribution Pend{z) of 
chain ends. We found that Pend{z) is very closely pro- 
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FIG. 7: (a) Rescaled values of the monomer density, [D + 
1) p{z) against ^ = z/(D + 1). Also plotted is the function 

/o(C) = 18.74 m-of'"'- 

(b) The data for D — 120 plotted against a modified scaling 
variable, = {z + S)/{D + 1 + 25), and divided by foi^s), for 
three different values of S. 

(c) The data for all D, again divided by fo{^s) with 5 — 0.15. 



portional to (^5(1 - ^s))°-^^^ with S = 0.3, but the (very 
small) deviations are highly significant. Taking them into 
account, we find 



(22) 



near the walls, with 8 
Eq. IHl). 



0.2. This agrees nicely with 
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FIG. 8: log-log plots of R^^{D)/N'^''^ versus D. The dashed 
line is 0.697 {D + l)-°-^^^. 



FIG. 10; Data collapse for testing the cross-over ansatz 
Eq. l2J for Domb-Joyce walks. 
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FIG. 9; Rescaled values of the probability pend{z) that the 
chain end is at the distance z from a wall, against (^s = {z + 
5)/{D -1-1-1-25) with 5 = 0.3. The solid line is the function 
ff(6)=4.78 (6(1 -6)^''. 



C. Domb-Joyce model 

Domb-Joyce chains with interaction strength w = 0.6, 
which is very close to critical streiiRth w* where leading 
corrections to scaling vanish [l^ Il3l |. were studied for 
slab widths up to D = 80. Chain lengths were up to 
N = 72000. The analysis of the data was done exactly 
as for the self avoiding walks described in Sec. IIIB. 

As expected, the scaling functions <&(??) and '5(77) are 
very similar to those for SAWs (see Figs. ^1 and 
This universality verifies that the amplitudes and criti- 
cal exponents discussed in Sec. II are essentially correct, 
although this should not be taken too serious: Such data 
collapse plots are not very sensitive to details (look at 
the huge range of scales in Fig. 1111 '). 

Estimates of the critical fugacities are shown in Fig. 1121 
where we plot 11 d — Moo against D + 1. The straight 
line, which again represents the extrapolation to large 
D, provides the estimate a = 0.2813(6). 

Plots of the monomer density profile (Fig. I13|l are very 
similar to those for SAWs. But the extrapolation length 
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FIG. 11: Data collapse for testing the cross-over ansatz 
Eq. for Domb-Joyce walks. 



is now much smaller, S « 0.04 as compared to 5 w 0.15 for 
SAWs. This is a first indication that corrections to scal- 
ing are indeed smaller in the DJ model. More important, 
also the scaling curve in Fig. 113b looks slightly different 
from that in Fig. [7|3: It is considerably smaller at the 
walls, with lim;,^o, d^oo D^+^/''^z-^/''^p{z)/A = 0.71(3) 
as compared to 0.87(4) for SAWs. Given the fact that 
scaling corrections should be smaller for the DJ model (in 
spite of the somewhat smaller values of D), we consider 
the DJ value as more correct, and blame the discrepancy 
onto scaling corrections for SAWs. With this new esti- 
mate of lim2_,o, D^oo Dp{z), and using the non- universal 
amplitude a determined earlier, we obtain our final esti- 
mate for the universal amplitude ratio B, 

B = 1.70 ± 0.08. (23) 

This is only 2 standard deviations away from the e- 
prediction B = 1.85 of Eisenriegler S], which we consider 
as good agreement. 

Finally, we do not show our data for Rn,\\ {D) /N'^^ and 
for the end monomer profile, since they are very similar to 
Figs. IHl and El But again the seemingly perfect agreement 
of the end point density profile is again as deceptive as 
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FIG. 12: Log-log plot oi hd ~ fJ,oo against D + 1 for the Domb- 
Joyce model. The dashed line is — Moo = 0.2813(1) + 
j^^-i/i's ^ith 1/3 used as constraint. 



FIG. 14: Rescaled values of the probability pend{z) that the 
chain end is at the distance z from a wall, divided by the 
function (7(^4) = 4.358(6(1-6)°-*', against 6 = {z+S)/{D+ 
1 + 25) with 5 = -0.02. 
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FIG. 13: (a) Rescaled values of the monomer density {D + 1 + 
2(5) p{z) of the Domb- Joyce model against £, = {z + S)/{D + 
1 + 25) with 5 = 0.06. Also plotted is the function fo{0 = 
18.74 (61 - O)'''"^ • (b) The same values as in panel (a), but 
divided by /o(6) with 5 = 0.04. 



it was for SAWs. This time our best estimate for the 
scaling of the end point distribution is 

Pendiz) ^ ^"-^KD , (24) 

with S « —0.02, in even better agreement with Eq. H18|) 
than the estimate for SAWs. The rescaled density pro- 



file divided by g{^s) = 4.358(^5(1 - 6))°'^^ is shown in 
Fig.Hl 



IV. SUMMARY 

We have presented high statistics simulations of 3-d 
polymers, modelled as walks on a simple cubic lattice 
with either hard or soft excluded volume interactions, 
with chain length up to 80000 on slabs of widths up to 
120. This was possible with the PERM algorithm with 
Markovian anticipation. The fact that PERM gives by 
default very precise estimates of free energies allowed us 
to measure precisely the forces exerted onto the walls, by 
measuring how the critical fugacities depend on the width 
of the slabs. We verified all critical scaling laws predicted 
for this problem, including the scaling of monomer and 
end point densities near the walls and the scaling of the 
total pressure with chain length and with slab width. 

The theoretical prediction most difficult to verify nu- 
merically concerns the amplitude ratio between the pres- 
sure onto the wall and the monomer density close to the 
wall. Previous simulations had not been able to obtain 
this with sufficient precision, and also in the present pa- 
per we had serious problems when using self avoiding 
walks with strict (hard-core) self repulsion. This might 
not be so surprising, given the well known fact that 
SAWs show rather large corrections to scaling. These 
corrections to scaling can be minimized by going over to 
Domb-Joyce polymers (characterized by soft repulsion) 
with carefully adjusted strength of the repulsion (sim- 
ilarly, for off-lattice bead-spring models, one can adjust 
the ratio between bead size and equilibrium spring length 
to minimize corrections to scaling) . It was only when go- 
ing over to this Domb-Joyce model, that we could verify 
in detail all theoretical predictions. 

Thus we have shown, first of all, that already the field 
theoretic e-expansion to first order in e, as implemented 
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in Ref. |3|, gives correct results. This was not obvious, 
in particular in view of persistent previous difficulties 
to verify it by Monte Carlo simulations. Secondly, we 
have demonstrated again the importance of using mod- 
els with minimized corrections to scaling. And last not 
least we have again shown that recursive sequential sam- 



pling methods with re-sampling (of which PERM is 
a particular implementation) can be very efficient. 
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